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Abstract 

We show that a wide variety of non-linear cellular automata (CAs) can 
be decomposed into a quasidirect product of linear ones. These CAs can 
be predicted by parallel circuits of depth 0(log 2 1 ) using gates with binary 
inputs, or 0(log t) depth if “sum mod p” gates with an unbounded number 
of inputs are allowed. Thus these CAs can be predicted by (idealized) 
parallel computers much faster than by explicit simulation, even though 
they are non-linear. 

This class includes any CA whose rule, when written as an algebra, is 
a solvable group. We also show that CAs based on nilpotent groups can 
be predicted in depth C>(logt) or 0(1) by circuits with binary or “sum 
mod p” gates respectively. 

We use these techniques to give an efficient algorithm for a CA rule 
which, like elementary CA rule 18, has diffusing defects that annihilate 
in pairs. This can be used to predict the motion of defects in rule 18 in 
C9(log 2 t) parallel time. 

PACS Keywords: 02.10, 02.70, 05.45, 46.10 


1 Introduction 

The direct product is a very basic notion in mathematics. Two algebras, dy¬ 
namical systems, or members of any other category can be paired so that their 
components act independently of each other: in algebras, (ai,6i) • (02,62) = 
(ai • 02,61 • 62). 
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More generally, we can have a quasidirect product where the first component 
is independent of the second, but not vice versa: 

(or, bi) • (<22,62) = (ai • (22, bi © ai> o 2 ^2) 

where ai and 02 determine some local operation ® ai ,a 2 011 the b’s. We will be 
particularly interested in the case where the b's form an Abelian group, and the 
local operations have the form 

^1 ©ai,a 2 b 2 fai,a 2 (bl') A <tai,a 2 (^2) A b (1) 

We will denote such a product by A® B. 

The semidirect product is a special case from group theory, in which g is the 
identity and / depends only on <22: 

(<2l, b±) • (<22, b 2 ) = ( a i ' °2, fa 2 (bi) ' b 2 ) 

For instance, if B is a normal subgroup of a group G (so that g~ 1 bg £ B for all 
g £ G and b £ B) and if A is a complement of B (so that every g £ G can be 
written uniquely as g = ab where a £ A and b £ B) then 

9 \ ‘ 92 = n\b\ ■ a 2 b 2 = <21(22 • {a^h a 2 )b 2 

and G is a semidirect product A® B with f a {b) = a~ 1 ba. 

This idea can be extended to dynamical systems. Rather than a direct 
product where two components evolve independently as in ( a',b') = $(a, b) = 
(. f(a),g(b )), we can have a quasidirect product of the form 

$(a, 6 ) = {f(a),g a (b)) 

The second component becomes a non-autonomous dynamical system, varying 
in time and space in a way controlled by the first component. If a dynamical 
system can be decomposed in this way, and if we have efficient algorithms to 
predict both / and the non-autonomous g, then we can predict the system as a 
whole. 

Cellular automata (CAs) are dynamical systems on the space of sequences 
over some finite alphabet, of the form 

^(<2)2 — • • • 1 a i, . • • , <2i-|-r) 

where r is the radius of the rule. By combining blocks of 2 r sites together, as 
shown in figure 1 , we can convert any CA into one with r' = 1/2 where each 
site has only two predecessors in a staggered space-time: 

$(a)j = 0(<2*—X/2, °2+ 1/2 ) 
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Figure 1: By blocking together 2 r sites, we can transform any CA into one on 
a staggered space-time with r' = 1/2. Here r = 2. 

(See |L9| for a study of this transformation’s algebraic properties.) We can then 
think of the CA rule as a binary algebra, where </>(a, b) = a ■ b or ab for short. 
The light-cone below an initial row becomes 


ao a i ci2 <23 

<2o<2l £ 11(22 02(23 

(a 0 a 1 )(a 1 a 2 ) (aia 2 )(a 2 a 3 ) 

((a 0 a 1 )(a 1 a 2 )) ((aia 2 )(a 2 a 3 )) 

and so on. Several authors have used this approach to explore CA properties 
such as partial reversibility |)J and periodicity Q. 

Predicting a cellular automaton t time-steps into the future is believed to 
be no easier in general than simulating it explicitly. To do this, we have to 
calculate all the CA states in a light-cone of depth t, which takes 0(t 2 ) serial 
computation steps (Q(t d+1 ) in d dimensions) or 0(t) in parallel. 

However, in Jfij we show that CAs whose rules satisfy various algebraic 
identities can be predicted in parallel time 0(log fe t) for some k , qualitatively 
faster than explicit simulation. We term these CAs quasi-linear: they are non¬ 
linear, but efficiently predictable nonetheless. 

It would be highly surprising if this were true for all CAs: since CAs exist 
which can simulate universal Turing machines (e.g. [fl-i| ), predicting a CA for a 
linear (or polynomial) amount of time is P-complete in general ||llf . Particular 
CA rules based on the Ising model or on majority-voting in three or more 
dimensions can also be shown to be P-complete [^8| . 

Therefore, it seems worthwhile to extend the class of quasi-linear CAs as far 
as possible, to explore what is probably a very rich hierarchy between linear 
dynamical systems and computationally universal ones. 

A preliminary version of these results appeared in [Q . 


2 Definitions 

An algebra (A, •) is a function from A x A to A, written a ■ b or simply ab. The 
order of an algebra is the number of elements in A. We will concern ourselves 
here with algebras of finite order. 
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The direct product A x B of two algebras is the set of pairs ( a,b ), with 
(oil bi)(a,2, 62) = (0102,6162). 

A quasigroup is an algebra whose multiplication table is a Latin square , in 
which every element of A occurs once in each row and each column. Quasigroups 
correspond to permutive CAs, in which 6) is a one-to-one function of each 
of its inputs (more generally, its leftmost and rightmost inputs) when the others 
are held fixed. 

An identity is an element 1 such that la = al = a for all a. An inverse of 
an element a is an element a -1 such that o _1 a = aa ^ 1 = 1. 

An algebra is associative if a(bc) = (ab)c for all a,b,c £ A. An associative 
algebra is called a semigroup. An associative quasigroup is a group. Groups 
have identities and inverses. 

An algebra is commutative if ab = ba for all a,b £ A. Commutative groups 
are called Abelian. The cyclic group Z p = { 0 , 1 ,... ,p — 1 } with addition mod p 
is Abelian. 

A function / on an algebra is a homomorphism if f(ab) = f(a)f(b). Homo- 
morphisms of Abelian groups can be represented as matrices. An isomorphism 
is a one-to-one and onto homomorphism; we will write A = B if A and B are 
isomorphic. An automorphism is an isomorphism from an algebra to itself. 

A subgroup B of a group A is a subset such that 6162 £ B for all 61,62 € B. 
The subgroup generated by all possible products of elements in a subset S is 
written ( S). 

A subgroup B is normal if a~ 1 ba £ B for all 6 £ B and all a £ A. For any 
normal subgroup B of A 1 there is a factor group A/B and a homomorphism 
from A to A/B that sends all elements of B to 1 . 

In a non-Abelian group, the commutator of a and 6 is [a, 6] = a~ 1 b~ 1 ab , so 
ab = ba[a, b]. The commutator subgroup G' = ([G, G]} of a group G is the set 
of all elements that can be written as products of commutators; then G/G' is 
Abelian. 

As a model of computation, we will consider families of circuits of varying 
depth and with different types of gates. The following classes of problems are 
those for which for all n, there is a circuit C n of depth O(\og k n ) and size 
polynomial in n that outputs the answer for inputs of size n: 

• NC fc if the gates are ANDs and ORs with binary inputs, 

• AC fe if the gates are ANDs and ORs with an unbounded number of inputs 
(unbounded fan-in), and 

• ACC fc [p] if the gates are ANDs, ORs, and “sum mod p” with an un¬ 
bounded number of inputs. ACC fe is the union U p ACC fe [p]. 

Then NC = UfcNC fc = UfcAC fc is the class of problems solvable in polylog- 
arithmic time by an idealized parallel computer with a polynomial number of 
processors, and is considered a good definition of problems that are efficiently 
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parallelizable. (A more realistic model of parallel computation would take the 
time and cost of communication between processors into account, assuming 
some finite-dimensional topology.) 

Families of circuits are uniform if there is a simple algorithm that generates 
C n when given n as input. Typically LOGSPACE-uniformity is used, i.e. a 
Turing machine that generates C n using d(logn) bits of memory. 

P is the class solvable by a deterministic Turing machine in polynomial time. 
It is easy to see that NCCP, but like P C NP this inclusion is believed, but 
not known, to be proper. From the definitions we have 


AC 0 c ACC 0 [ 2 ] c ACC 0 C NC 1 C AC 1 C ACC 1 C NC 2 C-CNCCP 

The parity function is clearly in ACC 0 [ 2 ], and Ajtai and Furst et al. 0 0 
have shown that it is not in AC 0 . Razboro v fl 25 |l has shown that majority is 
in NC 1 but not in ACC°[ 2 ], and Smolensky [f 28 | has shown that ACC°[p] and 
ACC°[g] are incomparable if p and q are distinct primes. Thus the first and 
second inclusions are proper. However, for all anyone has been able to prove, 
ACC 0 [6] could be equal to P (or even NP for that matter). 

A problem is P -complete if instances of any other problem in P can be 
converted to it by a LOGSPACE algorithm. It is generally believed that 
P-complete problems are inherently sequential , and cannot be efficiently paral¬ 
lelized; if any P-complete problem is in NC, then P = NC and all polynomial¬ 
time problems can be solved in polylogarithmic parallel time [jll], HU. 

The canonical P-complete problem is Circuit Value: what is the output 
of a given Boolean circuit, given the truth values of its inputs? Since truth 
values at each level of the circuit can affect those on the next level in arbitrary 
ways, it is hard to see how to get the output without going through the circuit 
level-by-level. 

An algebraic circuit over an algebra (A, •) is a circuit where each gate outputs 
the product a ■ b of its inputs, rather than implementing the standard Boolean 
functions; the Circuit Value problem for various classes of algebras has been 
studied by a number of authors 81 | 22 | . Predicting an r = 1/2 CA is clearly a 
special case of Circuit Value, where the circuit has a simple periodic structure 
in space and time. However, a number of the results we prove below for CAs 
will in fact hold for algebraic circuits of arbitrary shape. 

As shorthand, we will say a CA is in NC fc (AC fc , ACC fe ) if it can be 
predicted by circuits in these classes. Since a CA’s input consists of n = 2 rt+l = 
0 (t) initial sites, we will use n and t interchangeably in our O's. 


3 Non-autonomous additive CAs 

We will repeatedly use the fact that 

Lemma 0. The sum ofn elements of a finite Abelian group can be calculated 
in constant depth with ACC gates, or in 0 (logn) depth with binary gates. 
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Proof. Any finite Abelian group can be written as a direct product of 
cyclic groups G = Z P1 x Z P2 x • • ■ x Z Pfc . Addition of n elements of Z Pi can 
be carried out with a single ACC[pj] gate. Since a modpi = (aq mod Piq)/q, 
ACC[pi] gates can be simulated by ACC [pi.q] gates for any q, and in particular 
by ACC[p] gates where p = lcm(pi,p 2 , ■ • ■ ,Pk)- Alternately, a tree of binary 
gates of depth O(logn) can add the elements in pairs. H 

Now if a CA’s algebra is an Abelian group 4>(a , b) = a + b, a simple Green’s 
function method with Pascal’s triangle coefficients [jl7|, ^c| allows us to predict 
the CA in 0(t) serial time or 0(\ogt) parallel time, i.e. in NC 1 . In fact, this 
algorithm is in LOGSPACE-uniform ACC 0 since we can generate the f’th row 
of Pascal’s triangle using (D(\ogt) space. 

More generally, if a CA is of the form 

</>(ao, ffll, . . • , a2r) = /o®0 + Ml + • • ■ + /2r«2r + h, (2) 

where the ft are matrix-valued homomorphisms of an Abelian group ( A , +) and 
h is a constant element of A, we can represent the CA rule as a polynomial [|h| 

G(x) = fo + fix + f 2 x 2 H-b f 2r x 2r 

plus the constant h. Then we can think of the coefficients of G t (x) as the Pth 
row of a Green’s function for the CA, and for each t the final state is 

t t -1 t' 

s = y. G-ai+h y, y g\ 

2=0 t'= 0 2=0 

The multiplications can be carried out in constant depth, and by lemma 0 the 
sum is in ACC 0 . So any such CA is in ACC 0 . 

We now show that a non-autonomous version of (2) is still efficiently pre¬ 
dictable in parallel. 

Lemma 1. A non-autonomous CA of the form (2), where the ft are homo¬ 
morphisms of an Abelian group ( A , +) and h is an element of A, which vary in 
space and time independently of the CA state, is in ACC 1 . 

Proof. For simplicity, we will prove this for r = 1/2, where the CA rule is 
an operation of the form (1). A larger radius will simply increase the width of 
the light-cones, and the computation time, by a constant. 

Call the states in the light-cone St, x with t = 0 at the initial row, and s^o 
the leftmost state in the t’th row, as shown in figure 2. Then Sj iX ’s predecessors 
are St-\, x and and we write 

$t,x = — l,xi St— l,x-{-l) = ft,x(st— l,x) "b 9t,x(st — l,a;+l) “b 

(where we have replaced fi and f 2 with / and g for clarity) where ft tX and gt, x 
are homomorphisms of ( A , +) and is an element of A, all of which vary with 

t and x. 
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S 0,0 S 0,l S 0,2 


S 0,T 


S l,0 S l,l S 1,T-1 

S t-l.x S t-l,x+l 



S T,0 


Figure 2: The labelling scheme used in the text for the light-cone below the 
initial row. 


Now if ( tx') is in the light-cone below ( t , x), define C t ^ x '\t,x as the coefficient 
of St, x in St\ x '- Then we have 


„ f 1 if t' = t and x’ = x 

( ' X \ t,X \ ( ft',x ' 0 Ct' — l,x'\t,x) + (fft'.x' ° Ci'-l,x' + l|t,x) if t' > t 


Then it is straightforward to show inductively that the state st ,o at the bottom 
of a light-cone T steps high, with initial row sq.x for 0 < x < T, is 


T T T-t 

ST,0 = Ct,0|0,x(so,x) + EE Cr,o|t,x(/it,x) (3) 

x=0 t—1 x—0 

Since we are given the ht, Xl we can add all (T + 1)(T + 2)/2 of these terms 
together in depth O(logT) with binary gates or in constant depth with an 
ACC circuit by lemma 0; but we have to calculate the CT,o\t,x first. 

We will calculate the C t i <x '\t,x using a divide-and-conquer strategy. The 
influence of St, x on St\ x ' has to go through the intervening sites, so for any t" 
where t < t" < t 1 , we can write 

min (x,x' 

,x'\t,x 'y ^ ^t\x'\t'\x" ° ^t" ,x"\t,x (^) 

x // =max(x — t" +t ,x') 


as shown in figure 3. 

We can then use induction on increasing time intervals. Assume that the 
C t \ x '\t",x" and C t " tX "\t,x are known. Since the sum in (4) has at most (t'—t)/2+1 
terms and the compositions can be done simultaneously in constant time, by 
lemma 0 we can calculate C t f tX '\t,x i R depth 0(log(f' — t )) with binary gates or 
constant depth with ACC gates. 
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Figure 3: The divide-and-conquer strategy of equation 3. 


In particular, let t" = [{t' — t)/ 2J so that we double the time interval at 
each stage. Start by calculating (in parallel) all the C t ^ x f | t]X where t' — t = 1, 
then all those with t' — t = 2, and so on until we reach t' — t = T. This takes 
OlfogT) stages of induction, giving a total depth 0(}og 2 T) with binary gates 
or O(logT) with ACC gates, i.e. NC 2 or ACC 1 . 

Once we have CV.o|t,a; for all t,x we can easily add the sum in (3) as stated 
above to get the final state s-ro- So predicting the CA can be done in ACC 1 . 
■ 

We can generalize this further in two ways. First, for these sums to work, 
(A, +) simply needs to be commutative and associative, i.e. it can be a com¬ 
mutative semigroup rather than a group. Secondly, this algorithm works in any 
number of dimensions; the number of sites at t" between t and t! is proportional 
to ( t' — t) d , so the divide-and-conquer algorithm works just as fast. So in full 
generality we can state the following theorem: 

Theorem 1. Any complexity class of CAs containing ACC 1 is closed under 
quasidirect products with commutative semigroups where the local operations are 
of the form (1) where f and g are homomorphisms, or (2r + 1 )-ary operations 
of the form (2) where the fi are homomorphisms. Therefore, in any number of 
dimensions, suppose a CA is an iterated quasidirect product of the form 

C = (((C 0 ® Ci) ® C 2 ) • • •) ® C n 

where Co is in ACC 1 and the Ci for i > 0 are non-autonomous additive CAs 
of the form 


(fifai, 02,..., Ofc) — /i,i(ai) + /i j 2(«2) + • • • + /qfc(afc) + hi 



where the fij are homomorphisms of a commutative semigroup (Aj,+) and the 
hi are elements of Ai, depending on the state of Cj_i at each point in space-time. 
Then C is in ACC 1 . 

Proof. Since Co is in ACC 1 , we can calculate its state everywhere in the 
light-cone simultaneously with a circuit t 2 times as large. This and Lemma 1 
tells us how to calculate C\ everywhere in the light-cone, from which we can 
calculate C 2 everywhere, and so on. Iterate n times. H 

Example 1. The homomorphisms of Z 2 are zero and the identity, so the 
functions expressible as f{a) + g(b) + h are 0,1, a, a, b,b, a ® b and a © b. So 
any quasidirect product C ® Z 2 where the states of C select among these eight 
functions on Z 2 is in ACC 1 if C is. In particular, six quasigroups of order 4, 
including Z|, Z 4 , and four non-associative ones, can be written as quasidirect 
products Z 2 (g> Z 2 . 

Example 2. The automorphisms of Z 3 are /(a) = ±a. All quasigroups of 

3 elements can be expressed as ±a ± b + h. Therefore, any permutive CA which 
is a quasidirect product C® A, where A has three elements, is in ACC 1 if C is. 

Example 3. Suppose a CA has a vector-valued state (ai,..., a*,) at each 
site. Any r = 1/2 CA of the form 

(Pi(ai,b 1 ),P 2 (a 1 ,a 2 , h, b 2 ), ■■■, Pk(ai, ■ ■ ■, a k , b u ..., b k )) 

(or the obvious generalization for larger 7 ') is in ACC 1 if each Pi is linear in 
and bi but has arbitrary dependence on aj and bj for j < i. 

Since log(n c ) = O(logn). the divide-and-conquer algorithm in lemma 1 
actually reduces any polynomial size circuit of polynomial depth n c to an ACC 1 
circuit as long as the gates are functions of the form (2). Thus we can state the 
following: 

Corollary. The Circuit Value problem for algebraic circuits, over an 
algebra formed by an iterated quasidirect product of commutative semigroups 
where the local operations are of the form (2), is in ACC 1 . 

4 CAs based on solvable groups 

One interesting class of CAs consists of those for which <j>(a, b) = a ■ b is a 
non-Abelian group. Because of their non-commutativity these CAs are non¬ 
linear, i.e. they do not obey a principle of superposition, so Green’s function 
techniques don’t work. In [jnj we show an O(logt) algorithm for one such 
group, the Quaternions Qs , but other non-Abelian groups such as S 3 , the group 
of permutations of three elements, are left as open problems. 

We now show that a large class of finite groups have CAs in ACC 1 . First: 

Lemma 2. The set of algebras whose CAs can be predicted in a given amount 
of serial or parallel time (up to a multiplicative constant) is closed under finite 
direct products, subgroups, and homomorphisms. 
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Proof. For finite direct products G\ x G 2 x • • •, simply predict each of 
the Gi, either sequentially or in parallel. For subgroups, clearly an algorithm 
that predicts an algebra also predicts any of its subgroups. For a homomorphic 
image H = /(G), choose a pre-image f~ 1 (h ) for each element h £ H in the 
initial conditions, use the algorithm for G, and then apply / to return to H. 
Since / is a homomorphism, a simple induction shows that ..., h t ) = 

f{G\f-\h 0 ),...J-\h t ))) for all t. U 


In the language of universal algebra, this makes the set of CAs in a given 
complexity class a pseudovariety || . It would be wonderful if the pseudovariety 
corresponding to ACC 1 or NC, say, were finitely presented: that is, if a finite 
set of algebraic identities were necessary and sufficient for a CA to be in a given 
parallel complexity class. 

Now recall the following definition from group theory. 

Definition (12]. The derived series of a group G is the series of normal 
subgroups G = Go 3 Gi 3 ■ ■ • where Gj+i = ([Gi, Gi]) = G\ is the commutator 
subgroup of Gi. A group is solvable if its derived series ends in the identity {1} 
after a finite number of steps. Each of the factors Gj/Gj+i is Abelian. 

For instance, the derived series of S 4 , the group of permutations of 4 objects, 
is S 4 3 ^4 3 Zj 3 {1} where A 4 is the set of even permutations and consists 
of the identity and the permutations (12)(34), (13)(24) and (14)(23). 

Definition. A group is polyabelian if it can be written as a semidirect 
product of Abelian groups ((Ai<g)A 2 )(g)- • ')®A k with the parentheses associated 
to the left. (This could also be called an iterated split extension.) 

Clearly a CA based on a polyabelian group is in ACC 1 by theorem 1. Then: 

Theorem 2. Any CA whose algebra is a solvable group is in ACC 1 . 

Proof. We will show that any solvable group is isomorphic to a subgroup 
of a polyabelian group. Recall 129 that the wreath product AlB is a semidirect 
product B ® A B where A B is the set of functions a from B to A and elements 
of B permute their components. In other words, 


(b 1 ,a 1 )(b 2 ,a 2 ) = (hb 2 , f b2 (aq) oa 2 ) where f b {a)(x) = a(bx) 


The wreath product is useful for the following reason: if G has a normal sub¬ 
group N, then G is isomorphic to a subgroup of N l ( G/N ). 

Since every solvable group has an Abelian normal subgroup (the last non¬ 
trivial group in its derived series, whose commutator subgroup is { 1 }), by induc¬ 
tion it can be embedded in a wreath product of Abelian groups. In particular, 
if G = G 0 3 Gi 3 • • • 3 Gfc = {1}, then 

G C H k l {H k _ 1 \---{H 2 l Hfi)) 

where H k = G k -\/G k are the Abelian factor groups. This is a semidirect 
product of Abelian groups 

G C ((Ki ® K 2 ) • • • ® K k - 1) ® K k 
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where K\ = Hi and Ki +1 = H^\ Since G can be embedded in a polyabelian 
group, by lemma 2 a CA with G as its algebra is in ACC 1 . I 

As before, we can state the corollary for circuits in general: 

Corollary. The Circuit Value problem for algebraic circuits over solvable 
groups is in ACC 1 . 

For instance, the Quaternions Q s have the derived series Q s A {±1} A {1} 
since [a, b] = ±1 for any a,b and [±1,±1] = 1. The factor groups are Hi = 
Qs/{± 1} — ^2 and H 2 = {±1} — Z 2 , so Qs is a subgroup of Z 2 l Z 2 . This is a 

semidirect product Z| 0 Z| since zf 2 — Z 2 . 

In general, the Ki grow alarmingly in size, like • But most groups 

i times 

don’t require embedding in such large wreath products. Many small groups are 
themselves polyabelian, including the dihedral groups, groups of order p 3 for p 
an odd prime, any group of square-free order, all groups of order p 2 q where p and 
q are primes, and so on §@. All groups of order less than 32 are polyabelian 
except the dicyclic or generalized quaternion groups, which are factor groups of 
polyabelian groups twice their size, and the binary tetrahedral group of order 
24, which is a subgroup of a polyabelian group of order 12288 = 2 12 • 3 [|27]| . 

The smallest non-solvable group is A Fj . the simple group of order 60 (also 
called the icosahedral group). Since polyabelian groups are solvable (if G = 
((Ao 0 Ai) 0 A 2 ) 0 ..., then G' C [A\ 0 A 2 ) 0 ...) and since subgroups and 
factors of solvable groups are also solvable fiJJ, this group’s CA cannot be 
predicted by these methods. 

This leaves us with the following open question: is there an algorithm in 
ACC 1 , or NC fc for some k, for predicting CAs based on arbitrary finite groups? 
Barrington |j| and Beaudry et al. || have provided strong evidence to the 
contrary: since multiplication in A 5 can simulate NAND gates, circuits with 
non-solvable gates have P-complete Circuit Value problems. Unless there 
is some way to exploit the periodic structure of the circuit corresponding to a 
CA’s evolution, then, we would have 

Conjecture. CAs based on non-solvable groups such as A 5 are P- complete. 

Finally, we note that since the divide-and-conquer algorithm of lemma 1 ap¬ 
plies to circuits of arbitrary shape, we have actually shown that Circuit Value 
is in ACC 1 for solvable groups, or for polyabelian algebras formed inductively 
as in theorem 1. This is somewhat more general than the result of Beaudry et al. 
|j|, who show that Circuit Value is in the class DET for solvable semigroups, 
which is not known to be comparable with ACC 1 . 

5 CAs based on nilpotent groups 

We now show that a subset of the solvable groups have CAs in ACC 0 . 
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Definition Jl2| . The lower central series of a group G is the series of 
normal subgroups G = Ti 3 T 2 3 ••• where Tj+i = ([Tj, G]}. In other words, 
T 2 is the commutator subgroup G", T 3 is the subgroup generated by 3-element 
commutators [[a, 6 ],c], and so on. If Tfc + i = { 1 } for some k. we say that G is 
nilpotent of class k , and all commutators with more than k elements are 1. The 
nilpotent groups form a proper subset of the set of solvable groups. 

For instance, a nilpotent group of class 1 is simply an Abelian group. A 
nilpotent group of class 2 has commutators which commute with everything, so 
that [[a, 6 ],c] = 1 (such as the Quaternions where [a,b] = ±1 and [±l,c] = 1). 
And so on. 

Theorem 3. Any CA whose rule is a nilpotent group is in LOGSPACE- 
uniform ACC 0 . 

Proof. First, consider a nilpotent group of class 2. If its initial conditions 
are cio, a \,..., at, the leftmost two columns of its light-cone are 


a 0 a i 

a 0 ai aia 2 

ao afa 2 aia|a 3 

aoa\a\az ■ [ai, a 2 ] _1 01020304 • [o 2 , as ] -1 

The commutator arises since a 2 ai = Oia 2 [ai, a 2 ] , so 

( a 0 ala2)(aia 2 2 a 3 ) = aoa 3 a 2 [ai, a 2 ] _1 a2a 3 


Since commutators commute with everything we can move it to the right, leav¬ 
ing the ai in sorted order. Continuing in this way we get an “Abelian part” 

t (‘) 

rUo a x+i' times powers of commutators [ai,aj] where i < j. For t = 4,5,6 
these commutators are 

[oi, a 2 ] —4 [a 2 , o 3 ] —4 • [ai, a 3 ] — 1 

[a 1; a 2 ] _ 10 [a 2 , a 3 ] _ 24 [a 3 , a 4 ]~ 10 • [ai, a 3 ] _ 5 [a 2 , a 4 ]~ 5 • [ai,a 4 ] _1 

[ 01 , o 2 ] _ 2 °[a 2 , a 3 ] _ 84 [a 3 , a 4 ] _ 84 [o 4 , a 5 ]” 20 • 

[ai, a 3 ] _ 15 [a 2 , a 4 ] _ 35 [a 3 , a 5 ] -15 ■ [ai, a 4 ] _ 6 [o 2 , a 5 ] -6 • [a 4 , a 5 ] _1 


Since each site St, x in a light-cone is a product of St-i, x , which contains aj 


and st-i,x+i, which contains a}' * 1 , and since = afoj[aj, af] pq , we 
get (* Zl ) (ilZl 1) factors of I a . j . Oj] 1 in s t , x when we pass these powers of 
ai and aj through each other. Each of these commutators comes down to the 
final site st ,0 in {Z) wa y s > so 


st ,0 = I I a; 

0 <i<T 


n 


[aj, a,j 

0 <i<j<T 


_^T ^min(T-t,i —1) / T - t\ ( t - ( t- 1 \ 

I 2—j t = 0 i=max(0,j-t-l) ^ x Jyj — xjyi — x — 1) 
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(If [a, b] 2 = 1 for all a, b, this is just the algorithm for Q$ given in 0-) 

In general, if we have a nilpotent group of class k , then the final site can 
always be written 

-t,°=ir - n [[■ ■ ■ [ a u> ■ ■ •])(5) 

i i,j 

for some set of exponents c ,-(] 2 (1 < j < fc). This is clear by induction: 

by sorting the as* on the left, we generate 2-elenrent commutators [a,i,aj]. By 
sorting these we generate 3- and 4-element commutators, and so on, until we 
reach fc-element commutators which commute with everything. 

For each t, then, we have an expression (5) with 0(t k ) terms. Evaluating 
an expression of length m in a solvable group can be done in depth O(logm) 
with binary gates or constant depth with ACC gates |2). So nilpotent CAs are 
in ACC 0 . 

Finally, we show that this family of circuits is LOGSPACE-uniform, i.e. the 
circuit for predicting the CA t, time-steps into the future can be generated by an 
algorithm using only 0(log t) bits of memory. This is equivalent to calculating all 
the c^ i2 i .. Since j-element commutators are made by crossing commutators 

with fewer elements, each c ^ at St lX is a polynomial of ci’s at St- i, x and 
St-i, x +i where j' < j. These are polynomials of smaller commutators, and so 
on. 

So to get each , we have to calculate a tree of smaller ). This tree is of 
constant depth (at most k), and it ends in leaves of the form c,- 1 ^ = ^ which 

can be calculated in O(logt) space. So each of the can be calculated in 
logarithmic space, and since there are only polynomially many of them (which 
can be indexed with 0 (\ogt) bits) all of them can. I 

Since the same technique shows allows us to calculate the for any alge¬ 
braic circuit over a nilpotent group, we have the following corollary: 

Corollary. Any algebraic circuit over a nilpotent group is equivalent to an 
ACC 0 circuit. 

Note that this is not the same as saying that Circuit Value in general for 
circuits over nilpotent groups is in ACC 0 ; for this to be true, we would have to 
be able to take any circuit of arbitrary shape, and convert its topology, in con¬ 
stant parallel time, into the c^’s. We can say rather that any single circuit over 
a nilpotent group can be converted to an ACC 0 function of its inputs, and so 
e.g. LOGSPACE-uniform families of algebraic circuits become LOGSPACE- 
uniform families of ACC 0 functions. Then Circuit Value restricted to such 
a family is in ACC 0 , as CA prediction is. 


13 


6 Defining rules with words 

As another extension, suppose that we define a CA rule not simply as the 
product <f>(a 1 b) = ab in a given algebra A, but as the value of some word such 
as (f>(a,b ) = a 2 ba~ 1 . CAs with larger neighborhoods or in higher dimensions 
can also be defined this way, where 4>(a)i is the value of some word involving 
at and its neighbors as variables. We can also apply arbitrary homomorphisms 
to the variables in any way we please. Call such a CA rule word-defined on the 
original algebra A. 

Then since theorems 2 and 3 hold for circuits in general with solvable or 
nilpotent gates, and since each step of such a word-defined CA rule can be be 
written as a circuit with several gates of the original algebra, we can state the 
following: 

Corollary to Theorems 2 and 3. Any word-defined CA rule on a solvable 
or nilpotent group can be predicted in ACC 1 or ACC 0 respectively. 

Many interesting non-associative algebras can be written as words over non- 
Abelian groups; see [|| for a number of such constructions. 


7 Diffusing defect dynamics 

We close by showing that these techniques can in fact be applied to “classical” 
CA behaviors. Consider the following algebra on three symbols {0,1,1'}: 



0 

1 

T 

0 

0 

1 

T 

1 

1 

0 

0 

T 

T 

0 

0 


There are two subalgebras, {0,1} and {0,1'}, which appear as domains in the 
CA’s evolution: these are shown in black and gray in figure 4. Boundaries or 
defects between these domains act as diffusing particles, and annihilate when 
they meet in pairs. 

Theorem 4. The CA rule * shown in figure 4 is in ACC 1 . 

Proof. If we replace 0 with (0,0), 1 with (1,0), and 1' with (1,1), then * 
becomes a quasidirect product of T , 2 by Z2: 


(di, 61) * (<X2, 62) — (di © CZ2, f ai ,a 2 (h) © fiai,a 2 (^2) 

where / and g are homomorphisms depending on a\ and a 2 ' 

, th\ - i b if ai = 1, a 2 = 0 . _ J b if ai = 0, a 2 = 1 

J ai ,a 2 [t}) | 0 otherwise and 9a u a 2 [0) | Q otherwise 

Then by lemma 1 this CA is in ACC 1 . I 
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Figure 4: This CA rule has diffusing defects between domains (gray and black) 
that annihilate in pairs, similar to those in elementary CA rule 18. It can be 
predicted by ACC 1 circuits, and it matches rule 18 exactly for the motion of a 
single defect. 
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This rule is closely related to the two-state, r = 1 CA rule 18 (by the 
numbering system of |30|): 

Oi.iaiOi+i : 111 110 101 100 011 010 001 000 

a' : 00010010 

Blocking pairs of sites together gives the algebra 


• 

00 

01 

10 

11 

00 

00 

01 

10 

10 

01 

01 

00 

00 

00 

10 

10 

11 

00 

00 

11 

01 

00 

00 

00 


The subalgebras (00, 01} and {00,10} form two domains on which the CA acts 
like , and defects between them annihilate in pairs ill- If it weren’t for the 
product 10*01 = 11, the set (00, 01,10} would form a subalgebra isomorphic 
to * (with 00 = 0, 01 = 1, and 10 = 1') and 11 would not appear after the first 
time-step. 

However, if we choose initial conditions with a single defect, and choose the 
“phase” of the blocking rule so that we have (00, 01} on the left and {00,10} on 
the right, then the product 10 • 01 never occurs and the isomorphism is exact. 
We can use this approach to track multiple defects in parallel, as long as they 
don’t collide. Thus we have 

Corollary. The motion of a single defect in rule 18, or of multiple defects 
that do not collide, can be predicted in ACC 1 . 

By tracking single defects until they collide, we should be able to turn this 
into an algorithm for predicting rule 18 in polylogarithmic parallel time on 
average. However, there are rare worst cases, where a string of annihilations 
occur sequentially (each in the light-cone below the previous one) which appear 
difficult to parallelize. We leave this as an open problem. 

8 Conclusion 

Non-linear cellular automata that are easily predictable are akin to exactly solv¬ 
able models in statistical mechanics, or integrable non-linear partial differential 
equations. They help us categorize a hierarchy of dynamical systems between 
the purely linear and the computationally universal. One might hope that many 
systems in nature, such as turbulence, reaction-diffusion equations or flocking 
behavior might, at least in the abstract, be placed somewhere in this hierarchy. 

CAs based on groups, semigroups, quasigroups and so on are especially inter¬ 
esting, since they allow us to apply a rich and powerful theory of algebras to the 
problem of predicting the final state. We have shown here that CAs based on 
solvable and nilpotent groups can be predicted with ACC 1 and ACC° circuits 
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respectively, or in 0(log 2 t) and Oflogt) time by an idealized parallel computer 
whose gates have bounded fan-in. More generally, we have shown that the al¬ 
gebraic Circuit Value problem is in ACC 1 for solvable groups, and for a 
wider class of algebras formed by iterated quasidirect products of commutative 
semigroups. 

We have applied these results to a rule with diffusing defects, which can be 
used to predict the motion of defects in elementary CA rule 18. It would be very 
nice if some of the other complex nearest-neighbor CA rules, such as rule 22, 
30, or 110, could be analyzed this way; rule 30, for instance, has been proposed 
as a good cryptographic source of random numbers 131 , and a fast prediction 


algorithm for it would make this use inadvisable. Rule 55 has a more complex 
chemistry of travelling particles than rule 18 does Q, and it is unclear how to 
predict their motion efficiently. So far we have been unsuccessful at applying 
our methods to these CAs. 

Can these results be extended to CAs with looser algebraic structure? We 
note that of the 24 non-isomorphic quasigroups of order 4, for example, 14 are 
polyabelian and have CAs in ACC 0 . The other 10 are capable of expressing 
arbitrary Boolean functions and so their Circuit Value problem (and, we 
conjecture, their CA prediction problem) is P-complete & 

We hope these results add to existing P-completeness results for CAs such as 
majority-voting rules in three or more dimensions |li| as well as other discrete 
processes such as Ising dynamics and diffusion-limited aggregation [|jl| , and help 
flesh out this hierarchy from linearity to computational universality. 
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